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It is shown that a necessary and sufficient condition for an Archi- 
medean copula generator to generate a d-dimensional copula is that 
the generator is a d-monotone function. The class of d-dimensional 
Archimedean copulas is shown to coincide with the class of survival 
copulas of d- dimensional £i-norm symmetric distributions that place 
no point mass at the origin. The d-monotone Archimedean copula 
generators may be characterized using a little-known integral trans- 
form of WiUiamson [Duke Math. J. 23 (1956) 189-207] in an analo- 
gous manner to the well-known Bernstein-Widder characterization of 
completely monotone generators in terms of the Laplace transform. 
These insights allow the construction of new Archimedean copula 
families and provide a general solution to the problem of sampling 
multivariate Archimedean copulas. They also yield useful expressions 
for the d-dimensional Kendall function and Kendall's rank correlation 
coefficients and facilitate the derivation of results on the existence of 
densities and the description of singular components for Archimedean 
copulas. The existence of a sharp lower bound for Archimedean cop- 
ulas with respect to the positive lower orthant dependence ordering 
is shown. 

1. Introduction. Archimedean copulas are an important class of multi- 
variate dependence models which first appeared in the context of probabilis- 
tic metric spaces (see, e.g., [37]) and which enjoy considerable popularity in 
a number of practical applications. Following [22], any Archimedean copula 
has the simple algebraic form 

(1.1) C{ui, ...,Ud) = i^{^-\ui) + ■■■ + il^-\ud)), 

(«i,...,nrf)G[0,l]^ 
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where ip is a specific function known as a generator of C. As a consequence, 
many dependence properties of such copulas are relatively easy to establish 
because they reduce to analytical properties of the generator ip; see, for ex- 
ample, [11, 12, 15, 19, 26, 27]. Numerous statistically tractable parametric 
families of Archimedean copulas have been constructed and used in a vari- 
ety of practical applications, including multivariate survival analysis [5, 30], 
actuarial loss modeling [9, 21] and quantitative finance [4, 36]. 

While any d-dimensional Archimedean copula is necessarily of the form 
(1.1), the converse of this statement is not true. Somewhat surprisingly, 
the conditions under which a generator tp defines a d-dimensional copula by 
means of (1.1) have not yet been fully clarified except in two cases. Schweizer 
and Sklar [37] show that a generator ip induces a bivariate copula if and only 
if it is convex, whereas Kimberling [20] proves that ip defines an Archimedean 
copula in any dimension if and only if it is a completely monotone function, 
or equivalently, a Laplace transform of a nonnegative random variable. Kim- 
berling's condition is, however, not necessary for a given dimension d > 3 
and leads to limited dependence characteristics. Some authors [14, 26, 28] 
have required that tp only has derivatives up to order d which alternate in 
sign, a condition which, though sufficient and considerably weaker, is still 
not necessary. The present paper fills this gap by showing that the necessary 
and sufficient condition is that ip should have an analytical property known 
as d-monotonicity. This allows the possibility of d-dimensional Archimedean 
copulas without densities and reveals the existence of a sharp lower bound 
on the set of all d-dimensional Archimedean copulas with respect to the 
positive lower orthant dependence ordering as defined in [19]. 

The understanding of the necessary and sufficient conditions for ip leads 
to the insight that Archimedean copulas have a very natural geometric inter- 
pretation: they are the copulas that are implicit in the survival functions of 
d-dimensional ^i-norm symmetric distributions, also known as simplicially 
contoured distributions. These were first introduced by Fang and Fang [7] 
and comprise scale mixtures of the uniform distribution on the unit ^i-norm 
sphere. While it has been pointed out by Alfred Miiller that Archimede- 
an copulas can be constructed using ^i-norm symmetric distributions, this 
paper shows that these distributions are, in fact, absolutely central to the 
study of Archimedean copulas. 

Formalization of this link requires the consideration of a little-known in- 
tegral transform which appears in the work of Williamson [41] and which we 
call in this paper a Williamson d-transform. Essentially this transform plays 
the same role in the study of d-monotone generators that the Laplace trans- 
form plays in the study of completely monotone generators and knowledge 
of this link enables us both to propose a general solution to the problem 
of generating Archimedean copulas with arbitrary given generators and to 
construct and analyze a rich variety of new Archimedean copulas. 
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The paper is organized as follows: Section 2 discusses the necessary and 
sufficient conditions for a function tp to generate a d-dimensional Archi- 
medean copula. Section 3 establishes the connections between Archimedean 
copulas, £i-norm symmetric distributions and Williamson d-transforms. The 
theory is then utilized in Section 4 to prove several useful properties of Ar- 
chimedean copulas. This includes establishing conditions for the existence 
of densities and singular components, extending the notion of the Kendall 
function (a cornerstone of nonparametric inference for bivariate Archimede- 
an copulas introduced by Genest and Rivest [15]) to higher dimensions, and 
deriving a lower bound for Archimedean copulas with respect to the positive 
lower orthant dependence ordering. In Section 5 we comment on implications 
for stochastic simulation of Archimedean copulas and statistical inference. 

2. Archimedean copulas in a given dimension. This section gives nec- 
essary and sufficient conditions under which a generator ip induces an Ar- 
chimedean copula via (1.1). To begin, we need to establish basic notation. 
Throughout, x denotes a vector [xi, . . . , x^) in R"^; in particular, is the ori- 
gin. If not otherwise stated, all expressions such as x + y, max(x, y) or x < y 
are understood as componentwise operations. Furthermore, [x, y] refers to 
the set x ••• x [xd,yd] and M^|_ abbreviates the positive quadrant 

[0,00)''. Finally, ||x||i denotes the £i-norm of x, that is, ||x||i = J2i=i\xi\, 
and x+ denotes max(a;,0). 

The symbol X will be reserved for a random vector on M.'^ with distribution 
function H and survival function H, defined, respectively, by 

iJ(x) = P(X < x) and ^(x) = P(X > x) 

for any x € M'^. Note that H may be uniquely retrieved from H by means 
of the Sylvester-Poincare sieve formula and conversely. If the support of 
-ff is a subset of (0, 00)'^, H is uniquely given by its values in M^^ because 
obviously -ff (x) = ^(max(x, 0)), in which case we refrain from specifying the 
values of H outside of R^. Note, however, that if the associated probability 
distribution places mass on \ (0,(X))'^, in particular on the boundary of 
Wi^, H cannot be uniquely recovered from its restriction to M^J.. 

Finally, we will use difference operators defined as follows: Let / be an 
arbitrary d-place real function, x E M*^ and h > 0. Then the dth order dif- 
ference Ah/(x) is defined as 

Ah/(x) = Af^...Aij(x), 
where A^^, denotes the first-order difference operator given by 
^hj{^) = fi^U Xi-i,Xi + hi,Xi+i, ...,Xd) 

f i^l, . . . , Xi—l, Xi, , . . . , Xd) ■ 
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In the context of distribution functions, Ahi^(x) is the volume assigned by 
H to the interval (x, x + h] . This motivates us to define a function / : ^ ^ M, 
^ C M'^ fulfilling Ah/(x) > for any choice of x and h so that all vertices 
of (x, X + h] that lie in A as quasi-monotone on A. In other accounts quasi- 
monotone functions are referred to as d-increasing, but we prefer the former 
terminology here to avoid confusion with the notion of d-monotonicity of 
real functions, which will play a key role later on. 

Next, we recall few basic results which will be needed in subsequent dis- 
cussions. 

Definition 2.1. A (d-dimensional) copula is a function C:[0,1]'^^ 
[0, 1] satisfying 

(i) C(ui, . . . , Ud) = whenever Ui = for at least one i = 1, . . . , d. 

(ii) C{ui, . . . , Ud) = Ui if Uj = 1 for all J = 1, . . . , d and j / i. 

(iii) C is quasi-monotone on [0, 1]*^. 

Perhaps unfortunately, a survey of Archimedean copulas will amount to 
the investigation of certain probability distributions on given by their 
survival rather than their distribution functions. To avoid working with the 
more cumbersome formula for H in terms of it is convenient to formulate 
the following simple observation: 

Lemma 1. A d-place function H : — > [0, 1] is a survival function of a 
probability measure on M.'^ if and only if 

(i) H(—oo, . . . , — cxd) = 1 and H{x) =0 if Xi = oo for at least one i = 
l,...,d._ 

(ii) H is right- continuous, that is, for all x € M'^ it holds that 
Ve>0 3(5>0 Vy >x ||y - x||i < 5 ^ |^(y) - F(x)| < e. 

(iii) The function G given by G(x) = H{—j(.), x G W^, is quasi-monotone 
on W^. 

Proof. Since the proof is a standard exercise, we only provide a sketch. 
First, assume X is a random vector with survival function H. Then (i) and 
(ii) are immediate and (iii) is due to the fact that H{—x) = P(— X < x). Con- 
versely, suppose H satisfies conditions (i)-(iii). Then the right-continuous 
version G+ of G is a distribution function on R*^. Let X be a random vector 
with distribution function 6*+ and observe that H is the survival function 
of -X. □ 

Similarly, it will prove convenient to restate the original result by Sklar 
[38] in terms of survival functions. 
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Theorem 2.1. Let H be a d- dimensional survival function with margins 
Fi, i = 1, . . . ,d. Then there exists a copula C, referred to as the survival 
copula of H , such that 

(2.1) H{^) = C{h{xi),...,Fd[xd)) 

for any x € M'^. Furthermore, C is uniquely determined on D = {u G [0, 1]'^ : u G 
ran^i x ••• x ran^^} where ranFj denotes the range of Fi. In addition, for 
any u e 

C{M) = H{F{\ui),...,F^\ud)), 

where F~^ {ui) = m.i{x : Fi{x) < Ui} , i = l,...,d. Conversely, given a cop- 
ula C and univariate survival functions Fi, i = 1, . . . ,d, H defined by (2.1) 
is a d-dimensional survival function with marginals Fi, . . . ,Fii and survival 
copula C. 

In particular, if X is a random vector with survival function H and con- 
tinuous marginals Fi, . . . ,Fd and U is a random vector distributed as the 
survival copula C of H, we have that 

u4(Fi(Xi),...,Frf(Xd)) and X ^ {F^-^Ui), . . . , p-^Ud)). 

We are now in a position to turn our attention to Archimedean copulas. 
The latter were originally characterized by associativity and the property 
that C{u,u) < u for any u G [0, 1]; the present paper uses a more common 
definition based on the generator ip. 

Definition 2.2. A nonincreasing and continuous function ij:: [0,oo) — > 
[0,1] which satisfies the conditions V'(O) = 1 ^-iid \\mx^oo'il^{x) = and is 
strictly decreasing on [0, inf{x : -0(3;) = 0}) is called an Archimedean gen- 
erator. A d-dimensional copula C is called Archimedean if it permits the 
representation 

C(u) = V(V'"'(^i) + • • ■+i^~Hud)), u G [0, 1]<^ 

for some Archimedean generator ^ and its inverse : (0, 1] [0, oo) where, 
by convention, ■0(oo) = and (0) = inf : '>p{u) = 0}. 

Note that several authors define Archimedean copulas in terms of ■0"^ 
rather than ■0. The reason the above definition is favored here is that it 
leads, in this context, to simpler expressions, as will soon be apparent from 
the discussions below. 

A closer look at (1.1) readily reveals that 'ijj{ip~^{ui) + ■■■ -\- tp'^ (ud)) 
always satisfies the boundary conditions (i) and (ii) of Definition 2.1. There- 
fore, an Archimedean generator defines a d-dimensional copula via (1.1) 
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if and only if Tp{ip~^ (ui) + • • • + ip~^{ud)) is quasi-monotone. The following 
well-known result, which is Theorem 6.3.2 of [37], states conditions under 
which quasi-monotonicity holds in the bivariate case. 

Proposition 2.1. Let he an Archimedean generator in the sense of 
Definition 2.2. Then the function given by 

for u,v (z [0, 1] is a copula if and only if ip is convex. 

Proposition 2.1 does not extend to dimensions d > 3, as illustrated below. 

Example 2.1. It is a simple matter to check that the function ip{x) = 
max(l — x,0) is a convex Archimedean generator. However, 

ipii/j'^iui) H \- i/j'^ (ud)) = max(ui H \-Ud-d + 1,0), 

which is the Prechet-Hoeffding lower bound W{ui, . . . , Ud). As is well known, 
this function is not a copula for d>3, which can be seen by noting that it 
assigns negative mass to [1/2,1]*^. 

Consequently, stronger requirements on ip are needed. As will be shown in 
the sequel, these are based on the notion of d-monotone functions introduced 
below. 

Definition 2.3. A real function / is called d-monotone in (a, 6), where 
a, 6 G M and d>2, if it is differ entiable there up to the order d — 2 and the 
derivatives satisfy 

(-l)'=/«(x)>0, k = 0,l,...,d-2 

for any x G (a, 6) and further if (— l)'^"^ j(d-2) -g nonincreasing and convex 
in {a,b). Por d = 1, f is called 1-monotone in {a,h) if it is nonnegative 
and nonincreasing there. If / has derivatives of all orders in (a, h) and if 
>0 for any x in {a,b), then / is called completely monotone. 

Remark 2.1. Note that d-monotonicity as defined here is different from 
the notion of fc-monotonicity as used, for example, in approximation. More 
specifically, a function / : (a, 6) ^ M is called fc-monotone (or A;-convex in [33] 
and elsewhere) if its divided difference of order k is nonnegative. Without 
going into more detail at this point, we highlight the difference between the 
two concepts for k = 2: A;-monotone then means simply convex, while a 2- 
monotone function in our sense is nonnegative, nonincreasing and convex. In 
particular, therefore, a fc-monotone function is not necessarily /c'-monotone 
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for k' <k which contrasts with the above defined notion of d-monotonicity. 
Furthermore, it can be shown that, provided it exists, the /cth-order deriva- 
tive of a fc-monotone function is nonnegative, whereas in our case the deriva- 
tives alternate in sign. 

Definition 2.3 can be extended to functions on not necessarily open inter- 
vals. 

Definition 2.4. A real function / on an interval / C M is d-monotone 
(completely monotone) on /, d G N, if it is continuous there and if / re- 
stricted to the interior 1° of / is (i-monotone (completely monotone) on 

r. 

The central importance of d-monotonicity in the present paper is based on 
the following result, which relates d-monotonicity to the existence of survival 
functions. 

Proposition 2.2. Let f be a real function on [0, oo), p G [0,1] and H 
be specified by 

^(x) = /(II max(x, 0)||i) + (1 -p)l{x < 0}, x G M^. 

Then H is a survival function on M'* if and only if f is a d-monotone 
function on [0, oo) satisfying the boundary conditions lim^-^oo f{x) = and 
fiO)=p. 

The proof, together with that of Theorem 2.2 below, is deferred to Ap- 
pendix A. 

Theorem 2.2. Let tp be an Archimedean generator. Then C : [0, 1]*^ — > 
[0, 1] given by 

C(u) = V(V'"'(^i) + • • • + V'"'(^<i)), u G [0, 1]'^ 
is a d- dimensional copula if and only if tp is d-monotone on [0,oo). 

Note that Theorem 2.2 entails Proposition 2.1 as a special case. A further 
straightforward consequence is as follows: 

Corollary 2.1. Suppose tJj is an Archimedean generator which has 
derivatives up to order d on (0, cxd). Then tp generates an Archimedean cop- 
ula if and only if { — ^p^''\x) > for k = 1, . . . ,d. 

Obviously, d-monotonicity of an Archimedean generator ip does not imply 
that '0 is fc-monotone for k > d. In other words, d-monotone Archimedean 
generators do not necessarily generate Archimedean copulas in dimensions 
higher than d. This observation motivates the following notation: 
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Definition 2.5. denotes the class of d- monotone Archimedean gen- 
erators for d>2. In addition, ^'oo stands for the class of Archimedean gen- 
erators which can generate an Archimedean copula in any dimension d>2. 

On the other hand, if a function is d-monotone for some d>2, it is also 
/c-monotone for any 1 <k <d. This easy consequence of Definition 2.2 means 
that 

*2 5 *3 5 • • • • 

Before giving examples, which highlight in particular that \ ^d+i is 
nonempty, we note that the verification of the d-monotonicity of a generator 
tp can be considerably simplified by means of the following result, which is 
excerpted from [41], Theorem 4: 

Proposition 2.3. Let d>2 be an integer and ip be an Archimedean 
generator in the sense of Definition 2.2. Then ip is d-monotone on [0, oo) if 
and only exists on (0, oo) and is nonnegative, nonincreasing 

and convex there. 

For a complete account of the theory we can embed the well-known char- 
acterization of ^'oo due to [20] in our findings. 

Proposition 2.4. An Archimedean generator ip belongs to if md 
only if it is completely monotone on [0, oo). 

Example 2.2. Consider the generator -0^(3;) = (l — x)'^^ for some value 
of d>2. The L superscript anticipates a result in Section 4 in which this 
generator is shown to generate a lower bound for d-dimensional Archimedean 
copulas. We can verify that is d-monotone by computing {^]i)^'^~'^\x) = 
(— l)'^~^((i — 1)!(1 — x)+ and observing that (— 1)'^~^('(/'^)(''~^) is nonnegative, 
nonincreasing and convex on (0, oo). It is not, however, {d+ l)-monotone, 
since (^J')^'^~^^ does not exist at x = 1. Hence, V'rf ^ \ ^d+i- 

Example 2.3. The Clayton copula family has generator ipe{x) = (1 -|- 

-I In 

Ox)j^ . It is well known that for > Q this generator is completely mono- 
tone and can be used to construct a copula in any dimension. The case = 
should be understood as the limit (from either side) iPq{x) = Ximg^Qil + 

Ox)j^ = exp(— x), which generates the independence copula in any dimen- 
sion. 

The interesting case is < 0. When 9 = — l/(d — 1) for some integer d > 
2, then 'ipo{x) = ip^{9x) where ip^ is the generator in Example 2.2. The 
argument of that example holds and the multiplicative 9 is unimportant; 
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■061 generates copulas up to dimension d and, in fact, generates exactly the 
same copulas as ip^. For the general negative case let a = —1/9 > 0, write 

V'o(^) = i^-i/ai^) ^-rid observe first that il^a is convex or 2-monotone when 
a > 1. Moreover, when d>3 and a>d — l, then 

i-ir-'^^^^-'Hx) = a-('-'^ if (a - ^) (f - ^ 

fc=o ^ " 

exists and is nonnegative, nonincreasing and convex on (0,oo), whereas for 
other a values the derivative will either fail to exist everywhere or fail to have 
these properties; thus, by Proposition 2.3, the generator tpa is d-monotone if 
and only if a > d — 1. Summarizing, tpg is (i-monotone for a particular d>2 
if and only if 6* > - 1). 

3. Archimedean copulas and €i-norm symmetric distributions. For the 

remainder of the paper, assume that d>2. The beauty of Kimberling's char- 
acterization of ^'oo lies in its combination with the well-known Bernstein- 
Widder theorem (see, e.g., [40]). The latter states that an Archimedean 
generator is completely monotone on [0, oo) precisely when it is the Laplace 
transform of a nonnegative random variable, say W. In this case, as is well 
known, the d-dimensional Archimedean copula generated by is a survival 
copula of the survival function 

F(x) =V(||max(x,0)||i) = ^(e-ll'^^^(^'°)lli^) =^(e-^Sti'^^'^(^-0)), 

which is the survival function of the random vector X = y/W, where Y = 
(Yi, . . . ,Yd) is a vector of i.i.d. exponential variables, independent of W. X 
follows a mixed exponential distribution or frailty model of the kind that is 
often used to model dependent lifetimes; see [23, 29]. 

However, it is possible to give another representation for X. We can use 
the well-known facts that the random vector = Y/||Y||i has a uniform 
distribution on the d-dimensional simplex and and ||Y||i are independent 
to write X = RS^ where R = ||Y|[i/l^ and R is independent of S^. This 
shows that X follows a so-called ^i-norm symmetric distribution, a mixture 
of uniform distributions on simplices also known as a simplicially contoured 
distribution. 

Archimedean generators that are not completely monotone cannot appear 
as survival copulas of random vectors following frailty models, but they can 
all appear as survival copulas of random vectors following ^i-norm symmet- 
ric distributions. These distributions, introduced in [7], are formally defined 
as follows: 

Definition 3.1. A random vector X on M!^ = [0,00)^^ follows an £1- 
norm symmetric distribution if and only if there exists a nonnegative ran- 
dom variable R independent of where is a random vector distributed 
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uniformly on the unit simplex Sd, 



5d = {xGM^:||x||i = l}, 



so that X permits the stochastic representation 



X = RSd. 



The random variable R is referred to as the radial part of X and its distri- 
bution as the radial distribution. 

The general relationship between Archimedean copulas and £i-norm sym- 
metric distributions is quite similar to the relationship between Archimede- 
an copulas with completely monotone generators and frailty models outlined 
above. The difference is that the role of the Laplace transform is taken by 
another integral transform, which we call the Williamson d-transform and 
define as follows: 

Definition 3.2. Let X be a nonnegative random variable with distri- 
bution function F and d > 2 an integer. The Williamson d-transform of X 
is a real function on [0, oo) given by 



The class of functions which are Williamson d-transforms of nonnegative 
random variables will be denoted by W^. 

Remark 3.1. The Williamson d-transform belongs to the general class 
of Mellin-Stieltjes convolutions and is closely related to the Cesaro means. 
For more details, refer to [2], pages 194 and 246. 

It is clear that a Williamson d-transform always exists and is right-continuous 
at by the dominated convergence theorem. Far less obvious is the following 
result by Williamson [41]: 

Proposition 3.1. Let d>2 be an arbitrary integer. Then 

(i) yVd consists precisely of the real functions f on [0, oo) which are d- 
monotone on [0, oo) and satisfy the boundary conditions liuix^oo f{x) =0 as 
well as /(O) = p for p € [0, 1] . 

(ii) The distribution of a nonnegative random variable is uniquely given 
by its Williamson d-transform. If f = SU^F, then for x G [0, oo), F{x) = 
W^^f{x), where 



WdF{x)= f fl-yV 'dF{t) 

J{x,oq)\ tj 




if x>0 



if X = 0. 



2iJd V(x) = 1 - E 



d-2 



(-1) 



k=0 



kl 



id -1)1 
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Proof. The proof requires a reformulation of Williamson's original re- 
sult in terms of distribution functions. This is carried out in Appendix B. 

□ 

It may be noted that the proof of Proposition 3.1 indicates that if / € 
Wrf, then /(O) = p if and only if the corresponding distribution function 
F satisfies F{Q) = 1— p (i.e., has an atom at of size 1 —p when p < 1). 
Furthermore, since any / G Wd is strictly decreasing on [0, infjx : /(x) = 0}), 
we can conclude the following: 

Corollary 3.1. consists of the Williamson d-transforms of distri- 
bution functions F of nonnegative random variables satisfying F{0) =0. 

The role of the Williamson d-transform in the study of £i-norm symmetric 
distributions is indicated in Proposition 3.2 below, which also recalls two 
further properties of £i-norm symmetric distributions that will be used in 
this paper. These and many other results on £i-norm symmetric distributions 
have been derived by K. T. Fang, B. Q. Fang and collaborators; see the 
monograph [8]. 

Proposition 3.2. Let X = RSd be a random vector which follows an 
li-norm symmetric distribution with radial distribution function Fr. Then 

(i) The survival function H of^is given by 

(3.1) i?(x) = 2IJrfFR(||max(x,0)||i) + Ffi(0)l{x<0}, x e M'^. 

(ii) The density X exists if and only if R has a density. In that case, it 
is given by /i(|jx||i) = r(d) ||x||^~'^/ij(||x||i) where fn denotes the density of 
R. 

(iii) //P(X = O) = 0, thenR=\\X.\\i and = X/||X||i . 
Proof. See [8], Theorems 5.1, 5.4 and 5.5. □ 

The important implication of this proposition is that the survival func- 
tions of ^i-norm symmetric distributions reduce to 233d-Fij(||x||i), a simple 
function of ||x||i, whenever x € W^. In Proposition 3.3 below we extend this 
result to show the converse, namely that whenever a survival function is a 
function of ||x|ji on x G it must be the survival function of an £i-norm 
symmetric distribution, provided it places no mass on the boundary of M^^, 
except possibly at the origin. 

Proposition 3.3. Let be a random vector on M^|_. Then the following 
are equivalent: 
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(i) X has an ii-norm symmetric distribution. 

(ii) There exists a real function f on [0, co) so that the joint survival 
function H of ^ satisfies, for any x G M.'^, 

(3.2) ^(x) = /(II max(x, 0)||i) + (1 - /(0))l{x < 0}. 



Proof. The fact that an ^i-norm symmetric distribution satisfies (ii) is 
an immediate consequence of Proposition 3.2. The reversed imphcation can 
be estabhshed as fohows: First, note that if (ii) holds, / must be d-monotone 
on [0, oo) and satisfy /(O) S [0, 1] as well as lim2;^oo fix) = by Proposition 
2.2. Proposition 3.1 implies that / G Wd- Let be the distribution function 
Fr such that / = WdFR. Then H is precisely of the form (3.1) and X must 
have an £i-norm symmetric distribution with radial distribution Fpi. □ 



Remark 3.2. Proposition 3.3 is needed to fully establish the connection 
between £i-norm symmetric distributions and Archimedean copulas. In [7] 
a related result is proved that does not go quite far enough: the authors 
introduce a class Td of distributions on whose survival functions are of 
the form /(|[x|[i) for some function / on [0, oo) and argue that a member of 
Td is ^i-norm symmetric if / is d-monotone and /(O) = 1. What Proposition 
3.3 shows is that T^, in fact, coincides with the class of ^i-norm symmetric 
distributions as long as we simply exclude certain behavior on the boundary 
of M^. 

We are now in a position to bring all the pieces of the argument together 
in the following, which is the main result of this section: 

Theorem 3.1. (i) Let X have a d- dimensional ii-norm symmetric dis- 
tribution with radial distribution Fn satisfying Fii{0) = 0. Then X has an 
Archimedean survival copula with generator ij) = WdFjj . 

(ii) Let U be distributed according to the d-dimensional Archimedean 
copula C with generator tp. Then {tp~^ (Ui) , . . . , ip~^ (Ud)) has an ii-norm 
symmetric distribution with survival copula C and radial distribution Fr 
satisfying F^ = W^'^il), that is, 

^ ^-l)^xVW(x) {-ir-'x^-'i:t'\x) 

^'■'^ "^"^"^ = '"£0 (d^^- ' 

X G [0, 00). 

Proof, (i) From Proposition 3.2, part (i), it follows that the survival 
function of X satisfies -ff(x) = -0(11 max(x,0)||i) on M.'^ where -0 = SU^^-Fr. 
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By Corollary 3.1 we know that "0 ^ ^d- As shown in the proof of Theo- 
rem 2.2, since the marginal survival functions of H are the continuous func- 
tions Fi{x) = 'ip{ma,x{x,0)), we can conclude that the unique survival copula 
of X is the Archimedean copula generated by ip. 

(ii) The survival function of {^~^{Ui), . . . ,ip~^{U(i)) is given by H{xi, . . . , 
Xd) = C{'4){xi), . . . , ip{xd)) = V'dl max(x, 0)||i) on R'^. From Proposition 3.3 
we conclude that this is the survival function of an £i-norm symmetric dis- 
tribution and C is its copula. Since ip G ^'rf, it must be the Williamson 
d-transform of some distribution function of a nonnegative random vari- 
able R satisfying -Fr(O) = and this R is the radial part in the represen- 
tation {'^"^{Ui), Tp~\Ud)) = RSd. Finally, Fr = 2U^V follows by (ii) of 
Proposition 3.1. □ 

Remark 3.3. If we consider radial random variables with point mass at 
the origin and use these to construct £i-norm symmetric distributions, then 
Proposition 3.2, part (i), shows that these continue to have survival functions 
determined by the Williamson d-transform f{x) = WdFji{x), albeit in a case 
where /(O) < 1 and / ^ ^d- The ^i-norm symmetric distribution does not 
have a unique copula since it has discontinuous univariate margins. We enter 
the treacherous world of copulas for discrete random variables; see [13] . 

In light of the second part of Theorem 3.1 it will be useful to define the 
notion of an £i-norm symmetric distribution associated with a particular 
Archimedean generator. 

Definition 3.3. Let ip £ "^d and let Fr = SB^^V be the inverse 
Williamson d-transform of ^ specified by (3.3). Then Fr will be known 
as the radial distribution associated with ip in dimension d. If i? ~ -Fr is a 
random variable independent of S^, where is uniformly distributed on 5^, 
the distribution of RSd will be known as the ^i-norm symmetric distribution 
associated with ip in dimension d. 

Remark 3.4. Although there is a one-to-one relationship between ip and 
a particular radial distribution Fr and corresponding £i-norm symmetric 
distribution in a particular dimension, there is not a one-to-one relationship 
between the d-dimensional Archimedean copulas and the £i-norm symmetric 
distributions. Consider the radial variables R and R = kR for A: > 0. These 
have radial distribution functions FR{r) and Fj^{r) = FR{r/k), respectively, 
and give rise to different £i-norm symmetric distributions. Their Williamson 
d-transforms are ■0(a^) = ^dFR{x) and 'tp{x) = 'ip{x/k). However, these gen- 
erate exactly the same Archimedean copula in dimension d. 
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We close this section by giving examples of the use of Theorem 3.1 which 
show how we can calculate the radial distribution associated with a given 
Archimedean copula generator. 

Example 3.1. Consider the generator ipdi^) = (1 ~ ^)+~^ of the copula 
introduced in Example 2.2. Formula (3.3) easily yields that the radial dis- 
tribution associated with is, in fact, degenerate with R = l almost surely. 
Let Fj^ be the distribution function of another random variable R satisfying 
R = k almost surely, where k > 0. Theorem 3.1, part (i), implies that the d- 
dimensional ^i-norm symmetric distribution with radial distribution Fj^ has 
a survival copula with generator given by ip{x) = WdFj^{x) = (1 — x/k)'^^ . 
However, for any k > 0, this generates the same copula as ip^{x) = (1 — x)^^~^. 

Example 3.2. As is well known, the independence copula Il{ui, . . . , u^) 
is Archimedean with generator ip{x) = exp(— x). It is again easy to show 
by means of (3.3) that the radial distribution associated with this gener- 
ator is the Erlang distribution (or, alternatively, gamma distribution with 
parameter d) given by its density 

/R(x) = e-^-^-^— ^, xG[0,oo). 
The survival function of the associated £i-norm symmetric distribution is 

Ed 
y^^y i=i^'' by (3.1), which is obviously the survival function of in- 

dependent standard exponential variables. Note that, in conjunction with 
Proposition 3.2, part (iii), this shows that a simple way of generating a ran- 
dom vector Sd on the unit simplex is to generate d independent standard 
exponentials and divide each by the sum of the d variables. 

Example 3.3. Consider the Clayton copula family of Example 2.3 with 

-t /n 

generators ijjei^) = (1 + ■ Assume that a particular d > 2 and a 

particular 6 > —l/{d — l) are given. Suppose first that < 0, set a = —1/6* > 
and consider the generator V«(^) = i^-i/aix). Using (3.3) the associated 
radial distribution may be calculated to be 

d-l k-1 , rr\a-k 

(3-4) =l-E7fTn("-^-)(l-- 



1-E 

k=0 



a\ x\ I x 



k 



a J \ a 
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where (^) denotes the extended binomial coefficient. Moreover, Fji{x) = 
and Ffi{x) = 1 if, respectively, x <0 and x > a. Observe first that, if a = 
d—1, then this distribution simplifies to point mass at x = d — 1 as in Exam- 
ple 3.1, which is as we would expect. It is also an easy limiting calculation to 
establish that, as a ^ oo, the distribution function of R converges pointwise 
to the Erlang distribution in Example 3.2. 

If > 0, we can proceed in a similar manner by differentiating ipe{x) to 
show that for X G [0, oo) we have 

(3.5) Fn{x) = 1-1: nS^±:i!l^'=(i + fe)-(V^+fc). 

fc=o ^■ 

Again, as 6 ^0 this converges pointwise to the Erlang distribution function 
in Example 3.2. 

For illustration, consider the case when d = 3. li 6 = —0.5, the distribution 
of R is point mass at x = 2. For 6 > —0.5, R has a density and the copula 
is absolutely continuous. The radial densities are given in Figure 1 for the 
cases when 9 takes the values —0.3334, —0.3, —0.2, and 0.2. 



4. A new perspective on Archimedean copulas. The close relationship 
between Archimedean copulas and ^i-norm symmetric distributions delin- 
eated in Theorem 3.1 can be exploited in various ways. The important in- 
gredient is the fact that the distribution of the radial part of an ii- 
norm symmetric distribution can be explicitly retrieved from W^Fr and 
hence from H by means of (3.3), the inversion formula for Williamson d- 
transforms. From this point of view, working with Williamson d-transforms 
is more convenient than working with Laplace transforms for which the in- 
version formula cannot always be evaluated explicitly. 

The simple stochastic structure of £i-norm symmetric distributions can 
help us to explore analytical as well as dependence properties of Archime- 
dean copulas and several results of this kind are given in this section. 

4.1. Singular components of d- dimensional Archimedean copulas. The 
well-known Lebesgue decomposition theorem yields that any copula C can 
be written as 

C = Ca + Cs, 

where Ca and Cs are, respectively, the absolutely continuous and singular 
components of C. In other words, Ca is a distribution function of a (finite) 
measure on M'^ which is absolutely continuous w.r.t. the Lebesgue measure 
on M'^, that is, it has Lebesgue density. On the other hand, the measure 



16 



A. J. MCNEIL AND J. NESLEHOVA 




2 4 6 

X 



Fig. 1. Radial densities for the trivariate Clayton copula when 6 equals —0.3334, ~0.3, 
-0.2, and 0.2. 

induced by Cs is singular, that is, it is concentrated on a set of Lebesgue 
measure zero. 

Generally, the study of either component of a copula is not easy. However, 
in view of maximum likelihood estimation, it may be convenient to know 
whether a copula C has a density. In this section, we obtain more concrete 
results when C is Archimedean. Furthermore, singular components of Ar- 
chimedean copulas in the bivariate case are discussed in [12] and [27]; we 
extend these results to higher dimensions. 

To begin with, we should point out one potential fallacy. Let X denote 
a random vector with survival function H and continuous marginal sur- 
vival functions Fi, . . . ,Fd. Then, according to Sklar's theorem, there exists 
a unique survival copula C of H. Furthermore, as mentioned in Section 2, 
if U is a random vector whose distribution function is C, then 

\J = {Fi{Xi),...,Fd{Xd)) and X = {F^\Ui) , . . . , F^HUd)) . 

In other words, if and denote, respectively, the probability measures 
induced by C and H, then P*^ is an image measure of P^ with respect to 
a certain transformation and vice versa. Therefore, it may be tempting to 
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believe that if P is absolutely continuous, then so is P and vice versa. This 
may be true when Fi and F^^ obey the usual smoothness conditions required 
by the classical change of variable formula for Lebesgue integrals. In general, 
however, it is only guaranteed that Fi is differentiable almost everywhere and 
to investigate the relationship between the continuity properties of P^ and 

pH 

requires more involved change of variable formulas, which goes beyond 
the scope of this paper. Readers who are more interested in this subject 
should consult [17] for further information and references. We do provide a 
simple counterexample, however. 

Example 4.1. Take Fi^ to be the Cantor function (see [16], Example 
18.8). As is well known, Fn: [0,1] [0,1] is strictly increasing, continuous 
and satisfies F^(0) = as well as -^^(l) = 1 but F^(t) = a.e. in [0,1]. In 
particular, therefore, F^ provides an example of a singular univariate distri- 
bution function. Now, consider the following survival function H on M*^: 

d 

H{^) = '[[Fe{x,), xeM'^. 

4 = 1 

It is immediately clear that H does not induce an absolutely continuous 
probability measure. On the other hand, the marginals Ftf of H are con- 
tinuous and there exists a unique survival copula corresponding to H by 
Sklar's theorem. In fact, the latter is the independence copula 11 given by 
n(u) = nf=i ^i- Clearly, 11 has a density. 

In the case of Archimedean copulas, however, absolute continuity can be 
led back to absolute continuity of the associated ^i-norm symmetric distri- 
bution as follows: 

Proposition 4.1. Let C be a d- dimensional Archimedean copula with 
generator ^"^d- Let H stand for the distribution function of the ii-norm 
symmetric distribution associated with tp. Then 

(i) C is absolutely continuous if and only if H is. 

(ii) If tp G ^d+i; then C is absolutely continuous. 

Proof. Refer to Appendix C. □ 

It may be noted that Proposition 4.1 implies in particular that all lower- 
dimensional margins of an Archimedean copula are absolutely continuous. 

Proposition 3.2 states that an ^i-norm symmetric distribution is abso- 
lutely continuous if and only if its radial part is. Furthermore, the distri- 
bution function of the radial part is retrievable from ip using the inversion 
(3.3). These observations yield the following: 
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Proposition 4.2. Let C be a d- dimensional Archimedean copula with 
generator ip. Then C is absolutely continuous if and only if ip^"^'^^ exists 
and is absolutely continuous in (0,cx3). Furthermore, if C is absolutely con- 
tinuous, then its density c is given by 

^ ^W(V;-l(ui) + -.-+^-l(urf)) 

for almost all u S (0, 1)"'. 

Proof. Verification of this claim may again be found in Appendix C. 

□ 

Note that, in particular, the Archimedean generators considered in Corol- 
lary 2.1 obey the conditions of Proposition 4.2. 

Let us now summarize where we stand. If C is a d-dimensional Archime- 
dean copula, its generator ip belongs to either ^',^+1 or \ "^d+i- In the 
former case, C has a density by means of Proposition 4.1. In the latter, 
however, the density of C may or may not exist. Examples below illustrate 
situations that arise for ip (^^^X ^d+i- 

Example 4.2. tp^'^"^^ is absolutely continuous on (0, oo). In this case. 
Proposition 4.2 ensures that C has a density. An example in the bivariate 
case is as follows: Take tp on [0, oo) given by 

+ if xG [0,0.5), 

iP{x) = i ^a;2 - + f , if X G [0.5, 0.75), 
[o, if X G [0.75,cx)). 

Because ip'ix) = ^(x - 1) for x G [0,0.5) and ip'ix) = ^(2x - |) for x G 
[0.5,0.75), —ip' is nonincreasing and nonnegative on (0, cxo). In particular, 
therefore, ip €'^2- Furthermore, ip' is absolutely continuous in [0,00) and 
the bivariate Archimedean copula induced by ip is absolutely continuous. 
On the other hand, however, it can be easily verified that —ip' is not convex. 
Consequently, ip is not 3-monotone on (0, cxo) and hence not in ^3. 

Example 4.3. ip^'^~^'^ is continuous but not absolutely continuous on 
(0,00). We again provide an example in the bivariate case. Consider the 
function 

V.(x) = |i-^r^^(i-*)^*' ^^te[o,i], 

[0, ift>l. 
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where c = Jq F(r{l — t) dt and -F^ is the Cantor function introduced in Ex- 
ample 4.1. It is first clear that ip{0) = 1, V'(O) = and that ip is continuous 
on [0, cxo). Furthermore, 




if t > 1, 



for any x S [0, oo). In particular, —ip' is continuous, nonnegative and nonin- 
creasing and thus ip (z^2- However, ip' is singular on [0, 1]. It may be noted 
that the radial part of the associated ^i-norm symmetric distribution is also 
purely singular because, for almost all x € [0, 1], F'f,{x) = ^^^{^ — x) = 0. 

Example 4.4. tjj^' has jumps. Situations of this kind arise easily 
upon considering radial parts with atoms. To provide an example, set R 
to be a random variable independent of which follows the geometric 
distribution, that is, P(R = i) = p{l — for i € N. The random vector 

X = RSd is then £i-norm symmetrically distributed and its survival copula 
C is Archimedean with generator 

oo / \ d—l 

^P{x) = WdFR{x) = Y,{l--) p{l-py-\ xG[0,oo). 

In particular, because Fr has jumps and tp and ip^^^ k = 1, . . . ,d — 2 are 
continuous, (3.3) implies that ip^ has jumps and is thus not absolutely 
continuous. This confirms that C does not have density, which is already 
clear from the fact that R has atoms. 

So far, we have used the relationship between ^i-norm symmetric dis- 
tributions and Archimedean copulas in order to examine the existence of 
densities. As the rest of this section shows, the latter can be successfully 
employed for the investigation of the singular part of Archimedean copulas 
as well. We begin with a simple example: 

Example 4.5. A similar situation to Example 4.4 arises when we con- 
sider d = 2 and a radial part satisfying P{R = 1) = 2/3 and P{R = 2) = 1/3. 
The Archimedean generator of the survival copula of RS2 is then 

V'(x) = mFRix) = ^(1 -x)+ + ^(^i- 0^. 

By the same argument as in Example 4.4, the bivariate Archimedean cop- 
ula C induced by tp does not have a density. Figure 2 illustrates that the 
^i-norm symmetric distribution with radial part R is concentrated on two 
simplices, 

52(1) = {xGM^: ||x||i = 1} and ^2(2) = {x G : |[x|[i = 2}. 
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Fig. 2. 10,000 simulated values from the bivariate li-norm symmetric distribution of 
Example 4-5 (left panel). The middle panel shows the corresponding survival copula and 
the right panel shows its generator ijj. 



C is also purely singular as can be seen in the middle panel of Figure 2. It 
is not difficult to verify that the support of C is U A2 , where 

Ai = {{ui,U2)e[0,lf:ip-\ui)+tlr\u2) = i}, i = l,2. 

Note that if U ~ C, then P(U G ^1) = 2/3 and P(U G A2) = 1/3. 

The findings of Example 4.5 can be generalized. To do so, we first need 
to introduce additional notation. A level set L{s) of a copula C is given by 

L(s) = {uG [0,l]^:C(u) = s} 

for s G [0, 1]. For a d-dimensional Archimedean copula the level sets take the 
form 



L{s) 



{u,) = ^-\s)\, ifsG(0,l], 



uG[0,l]^:5]V'-i(t 
I 1=1 

|u G [0, 1]^ : ^ V^"i (n,) > (0) | , if s = 0. 



Proposition 4.3 below, which is a high-dimensional extension of the results 
due to [12] and Alsina, Frank and Schweizer (see [27], Section 4.3), indicates 
the mass placed by an Archimedean copula on its level sets. 

Proposition 4.3. Let C be a d- dimensional Archimedean copula with 
generator and let stand for the probability measure on [0, l]'^ induced 
by C . Then 



P^(L(.)) 



(-iOVlM^^/,('^-i) 



{d-l)\ 
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for s G (0, 1]. Furthermore, if ip ^(0) = oo, then P (L(0)) = 0. Otherwise, 

Pcim) = P^ju e [0, 1]'^ : X: ^~\ui) = V'-^o)! 

Proof. Refer to Appendix C. □ 

4.2. The Kendall function in d dimensions. For a bivariate Archime- 
dean copula C and a random vector {Ui,U2) = C, the Kendall function 
-f^c (also known as the bivariate probability integral transform) defined as 
the distribution function of C{Ui,U2) is a cornerstone of nonparametric 
inference for Archimedean copulas as introduced in [15]. Although this paper 
is not devoted to estimation, it is nonetheless worth noting that Theorem 3.1 
allows an easy generalization of several of the findings about Kq to higher 
dimensions. Results from this section also nicely complement those in [1], 
Example 3. 

Let C be a d-dimensional Archimedean copula with generator ip and U 
C be a random vector. The function Kq is, in analogy to the bivariate case, 
given as the distribution function of the random variable C(U), that is, 

Kcix) = P{C{Ui,...,Ud)<x). 

Clearly, C(U) is concentrated on [0,1]. For bivariate Archimedean copulas, 
Kc{x) can be given explicitly in terms of the generator and its left-hand 
derivative (see [15], Proposition 1.1). To establish an analogous result in the 
d-dimensional case, observe first that 

P(C([/i, ...,Ud)<x) = p(Y\^ i;~\U,) > ^l;~Hx)^ ■ 

Furthermore, let X be a random vector which follows the £i-norm symmetric 
distribution associated with ip. Sklar's theorem then implies 

xi(^"i(c/i),...,V-^(f/,)). 

In particular. Proposition 3.2 yields that the radial part P of X satisfies 

R^\\ni=T.^-'m 

1=1 

and is independent of X/|[X||i. The consequences are now immediate. 
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Proposition 4.4. Let C be a d- dimensional Archimedean copula with 
generator ip and let JJ be a random vector with U ~ C. Then 



d 



J2^~HU^) and 



i=l 



are independent. Moreover, for any j = 1, . . . ,d, the random variable 



is standard uniform. 



Proof. What remains to be shown is the uniformity of Vj , j = 1, . . . ,d. 
This follows easily by observing the fact that the univariate marginals of the 
uniform distribution on the simplex Sd follow the Beta(l, d— 1) distribution, 
that is, their survival functions are given by il^^{x) = (1 — x)"^^. □ 

Proposition 4.5. Let C be a d-dimensional Archimedean copula gen- 
erated by ip. Then 

r (-i)^-.(,-(o))^- ^,_., 
[d-iy. 

Kc{x) = { 2^ 

k=0 

, {-lr~\^^-\x)r~'^^^-'\^p~\x)) 

Proof. Because Kc{x) = P{R > ip~^{x)) where R denotes the radial 
part of the ^i-norm symmetric distribution associated with ip, the assertion 
follows directly from part (ii) of Theorem 3.1 and from Proposition 4.3. □ 

4.3. Archimedean copulas are bounded below with respect to the PLOD 
ordering. Recall that the positive lower orthant dependence (PLOD) or- 
dering is a partial ordering on the set of all d-dimensional copulas given as 
follows [19]: 

Ci -<cL C2 ^ Ci{u) < C2{u) for any u G [0, 1]'^. 

It is a well-established fact that any bivariate copula C satisfies W C 
where W is the Frechet-Hoeffding lower bound specified in Example 2.1. In 
dimensions greater than two, it still holds that W -<cL C, with the unfortu- 
nate exception that W is no longer a copula. Furthermore, it can be shown 
that for d > 3, no sharp lower bound of the set of all copulas with respect 
to the PLOD ordering exists (see [27], Theorem 2.10.13). 
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The situation for Archimedean copulas is different, however. As is well 
known, an Archimedean copula C whose generator is completely monotone 
is positive lower orthant dependent, that is, 11 C where 11 is the indepen- 
dence copula given by n(u) = nf=i (see [27], Corollary 4.6.3.). Proposition 
4.6 below shows that a similar statement can also be made in the general 
case. 

Proposition 4.6. Let C he a d- dimensional Archimedean copula. Then 

Cd -<cL C, 

where is the d-dimensional Archimedean copula with generator ip^ of 
Example 2.2. 

Proof. Refer to Appendix C. □ 

Figure 3 below illustrates the lower bound in dimension d = 3. 

A useful implication of Proposition 4.6 is that if C is the bivariate marginal 
distribution function of a d-dimensional Archimedean copula with generator 
ip (z^d and p an arbitrary measure of concordance in the sense of [35] , then 

p{Cy)<p{C), 

where p{C^) refers to p applied to the bivariate marginal distribution of C^. 
This applies in particular to Kendall's r, which is a concordance measure 
given by 

r(C)=4 [' [' Ciui,U2)dCiuuU2) - 1 = 4EiC{Ui,U2)) - 1, 
Jo Jo 

where {Ui,U2) refers to a random vector distributed according to C. 




Fig. 3. The left panel shows the support of the lower bound Archimedean copula Cj" for 
d — 3. Two thousand simulated points from the bivariate margin of are depicted in the 
right panel. 
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For bivariate Archimedean copulas, it was discovered in [11] that t(C) 
can be expressed in terms of the corresponding Archimedean generator, 



where the second expression results from the first by a change of variable. 
The following simple observation offers a geometric interpretation in terms 
of the associated £i-norm symmetric distribution: 

Proposition 4.7. Let C be a bivariate Archimedean copula with gen- 
erator ip and R be the radial part of the ii-norm symmetric distribution 
corresponding to C in the sense of Definition 3.3. Then 



Proof. Let (f/i, U2) be a random vector distributed as C. The assertion 

is then an immediate consequence of the fact that ^(■(/'"^(C^i) +^~^(C/2)) = 
tp{R), which has been established in Section 4.2. □ 

The lower bound on Kendall's tau t(C) is now straightforward. 

Corollary 4.1. Let C be the bivariate margin of a d-dimensional Ar- 
chimedean copula with generator ip. Then t{C) > —l/{2d — 3). 

Proof. The proof follows from (4.1) using elementary calculus. □ 

5. Practical implications. The results and connections established in this 
paper have a number of practical implications for working with Archimedean 
copulas. We confine ourselves to brief comments on implications for copula 
construction, stochastic simulation and statistical inference. 

5.1. Construction of new distributions. Many families of bivariate Archi- 
medean copulas have been proposed, but models that extend to dimensions 
three and higher are much less common. The majority of such copula fam- 
ilies have generators which are Laplace transforms of nonnegative random 
variables and, as discussed in Section 4.3, can only capture positive depen- 
dence. Examples of multivariate Archimedean copulas whose generators are 
not completely monotone seem to be rare; the standard examples are the 
Frank or the Clayton family introduced in Example 2.3. Theorem 3.1 offers 
an elegant solution to this problem: Any d-dimensional Archimedean cop- 
ula can be obtained simply as a survival copula of a d-dimensional £i-norm 



(4.1) 





r(C7) = 4S(^(i?)) - 1. 
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symmetric distribution with some radial part R. Concretely, by selecting 
a p-parametric family Tg = {Fg : ^ G C RP} of distributions on [0, oo), we 
obtain a parametric family Cg of d-dimensional Archimedean copulas whose 
generators are given by 

{M^)=^dFgix):eee}. 



Example 5.1. Consider the two parametric family of distributions with 
densities 



fa,b{x) 




if a < X < 6, 
otherwise, 

where < a < 6. fa,b{x) is the density of the reciprocal of a uniform distribu- 
tion on the interval a~^]. The survival copula of the ^i-norm symmetric 
distribution with this radial distribution has generator 

^ ah of xY~^ , 
1 - - dt 



b — a V t 



ab ((^ x\ / X 



xd{b — a)\\ b J ^ \ a 

Clearly, tpa,b{x) is not completely monotone because it does not have deriva- 
tives of all orders. In Figure 4 we show four examples of 2000 points gener- 
ated from the corresponding bivariate copula for a = 1 and different values 
of b. 



5.2. Stochastic simulation. For a given Archimedean generator ijj G 
generating a copula C in dimension d, it is of interest to have a general 
algorithm for generating random vectors U with distribution function C. If 
ip € ^'oo) then a known method of generation is to use the mixed exponential 
or frailty representation discussed at the beginning of Section 3. This idea 
can be traced back to [23]; see also [24]. A key feature of the method is that 
we need to invert the Laplace transform -0 to find a probability distribution 
function Fy/ on the positive half-axis, from which we must then be able to 
sample; this is not always straightforward. 

The insights of Section 3 give an alternative method. The steps of this 
method may also be found in [39], but are only justified there for generators 
■0 £ ^oo- In fact, the method may in principle be applied to any generator 
ip G for d>2. To generate a random vector from the d-dimensional copula 
C with generator ip, the algorithm is as follows: 
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a = 1 ; b = 1 .2 



a = 1; b = 2 




0,0 



0.2 



0.4 0.6 
U1 

a= 1; b = 5 



0.8 



1.0 



0.0 



0.2 



0.4 0.5 



o.e 



1.0 




0.0 0.2 0.4 0.6 0.8 
U1 

a = 1; b = 1000 



1.0 




r ■--•.■-*4,.-;- -j.'.:7i 





0.0 



o.z 



0.4 0.6 
Ul 



0.8 



1.0 



Fig. 4. Four examples of 2000 points generated from the copula of Example 5. 1 for a - 
and different values of b. 



1. Generate a random vector S = (^i, . . . , Sd) uniformly distributed on the 
d-dimensional simplex Sd- Based on Example 3.2, we can generate a vec- 
tor of i.i.d. standard exponential variates Y = (Yi,...,!^) and return 
5, = y,/||Y||i for i = l,...,(i. 

2. Generate a univariate random variable R having the radial distribution 

associated with ^ in the sense of Definition 3.3; this will generally 
require us to compute the inverse Williamson d-transform Fji = 'W^^iIj. 

3. Return U = (C/i, . . . , Ud) where Ui = ip{RSi) for i = 1, . . . , d. 
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The second step is the practically challenging part of the algorithm. In 
general we would use the inversion method to generate R although this may 
require both accurate evaluation of -Fr = SU^^V and numerical inversion 
of Fpt', see Example 3.3 for the distribution function that would have to 
be computed and inverted for the Clayton copula. An example where the 
method works particularly efficiently and allows us to generate samples from 
an attractive bivariate family is given below. 

Example 5.2. Consider the generator ip{t) = (1 — t^^^)+ where 9 >1; 
this is the second family considered in [27] (see Table 4.1). Since ip is convex, 
it follows that £ ^2 but tp ^ ^3 because ip^^^t) is not defined at t = 1. 

Since ip^^\t) = -e~^t^/'^~^ for t G [0,1) and V'+^t) = otherwise, we can 
use Theorem 3.1, part (ii), to show that 

U, ifx>l, 

which has a jump at x = 1 and yields a bivariate copula with a singular 
component. The distribution of R is easily sampled by inversion. In Fig- 
ure 5 we show four examples of 2000 points generated from this copula for 
different values of 6. Obviously, 9 = 1 corresponds to the Frechet-Hoeffding 
lower bound copula generated by ip^ and, as 9 ^ 00, we approach perfect 
positive dependence. It is easy to use Proposition 4.7 to calculate that the 
Kendall's rank correlation of this copula is given by r(C) = 1 — 2/9, which 
implies that the points in the second picture (9 = 2) are taken from a copula 
with a Kendall's rank correlation equal to zero. 

5.3. Statistical inference. Proposition 4.4 can be used to devise diagnos- 
tic tests for particular hypothesized copulas. Suppose we have d-dimensional 
data vectors Ui , . . . , that are believed to form a random sample from 
an Archimedean copula with generator ■0. (Alternatively they might be ob- 
servations from the empirical copula of some multivariate data, derived by 
applying the marginal empirical distribution functions to the data as in [10].) 
For each vector observation = (Uki, ■ ■ ■ , Ukd) we can form the quantities 
Rk = Eti i^-HU^) and Sfc = ii^-HUki), . . .,^p-HUkd))/Rk- 

The Rk data should be compatible with the radial distribution = 
W^^ip, whereas the data should be compatible with the hypotheses of 
uniformity on iS,^ and independence from the radial parts. Clearly, there 
are various possible ways of deriving numerical or graphical tests of these 
hypotheses. 

It is also possible to reduce the estimation of Archimedean copulas to 
a one-dimensional problem by focusing on radial distributions. Ideas devel- 
oped by [15] for the bivariate case may be extended to the general multivari- 
ate case. In that paper a nonparametric margin-free estimation procedure 
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0.2 0.4 0.6 0.8 1,0 0.0 2 0.4 0.6 0.8 1.0 

Fig. 5. Four examples of 2000 points generated from the copula of Example 5.2 for 
different values of 9. 

is proposed for the distribution function Kc{x) = P{C{Ui, U2) < x), where 
{1/1,112) are random variables distributed according to the unknown copula 
C. This procedure may be applied in the higher-dimensional case to de- 
rive a nonparametric estimate of Kc{x) = P{C{Ui, . . . , Ud) < x), which is 
also given by Kc{x) = P{R > tlj~^{x)), where R is the radial part of the 
^i-norm symmetric distribution associated with the generator ip of C. Thus 
the distribution Kq may also be computed under parametric assumptions 
for the copula generator ^ using Proposition 4.5. This makes it possible to 



MULTIVARIATE ARCHIMEDEAN COPULAS 



29 



estimate the copula by choosing the parameter values that give the best cor- 
respondence between the parametric model for Kq and the nonparametric 
estimate. 

APPENDIX A: PROOFS FROM SECTION 2 

Theorem 2.2 relies on several results on c?-monotone functions. First is 
the following lemma which summarizes Theorems 4, 5 and 6 of Chapter IV 
of [40]. 

Lemma 2. Let d>l he an integer and f a nonnegative real function on 
{a,b), a,6sM. // / satisfies 

(A.l) (A„)V(^)>0 

for any k = 1, . . . ,d, any x £ (a, b) and any h > so that x + kh £ (a, b), 
then 

(i) / is nondecreasing on (a, 6). 

(ii) If d>2, f is convex and continuous on {a,b). 

(iii) If d>2, the left-hand and right-hand derivatives of f exist every- 
where in (a, b). Moreover, for any a < x < y < b, 

fUx)<f'+{x)<fUy)- 

(iv) For any k = 1, . . . ,d—l and any a < x < y < b, {Afi)^f[x) < {Ah)^f[y) 
whenever h is chosen small enough so that y -\-kh < b. 

The second result which has proved of key importance is the following 
characterization of d-monotone functions: 

Proposition A.l. Let f he a real function on {a,b), a,beR and let f 
denote a function on (—6, —a) given hy f{x) = f{—x). Further, let d>l he 
an integer. Then the following statements are equivalent: 

(i) / is d-monotone in {a,b). 

(ii) / is nonnegative and satisfies, for any k = 1, . . . ,d, any x G (—6, —a) 
and any hi> 0, i = 1, . . . ,k so that (x + /ii + • • • + hk) € {—b, —a), 

Ah,---Ahj{x)>0, 

where A/^^. • • • A^j^ denotes a sequential application of the first-order differ- 
ence operator Ah given hy Ahf{x) = f{x + /i) — f{x) whenever x,x -\- h 
{-b-a). 

(iii) / is nonnegative and satisfies, for any k = 1, . . . ,d, any x £ (—6, —a) 
and any h> so that x -\- kh £ (—6, —a), 

{Ah)'f{x)>0, 

where {A^)^ denotes the k-monotone sequential use of the operator A^. 
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Proof. 

Proof of (ii) (iii). First note that (iii) is a special case of (ii). For the 
reverse imphcation, assume without the loss of generality that d>2 and fix 
an arbitrary x in (—6, —a) and hi > 0, i = 1, . . . ,d such that x + hi + ■ ■ ■ + ha G 
(—6, —a). Now, let fi, £ = 1, . . . ,d, denote a function on (—6, —a — {hd + • • • + 
hn)) given by 

h{y) = Afc^ • • • ^hj{y), y G {-h, -a-{hd + --- + he)). 

Observe in particular that any fe, i = 1, . . . ,d, is nonnegative. One can even 
show that the functions fe, £ = 2, . . . ,d, satisfy (A.l) of Lemma 2 for any 
k = l,...,£-l. 

This claim can be verified by induction: because / satisfies the assump- 
tions of Lemma 2, the assertion (iv) thereof yields that (A/j)^'/(y) < {Ah)^f{y + 
hd) for any k = 1,. . . ,d — 1 whenever h is small enough so that y + kh ^ 
{—b, —a — hd)- But 

{A,)'^f{y + hd) - (A,,) V~(y) = (A,,) Vd(y) 

for any k = 1, . . . ,d — 1 and hence the beginning of the induction is estab- 
lished. 

Now suppose that satisfies (A.l) for any k = 1, . . . ,£. The fact that 
My) = fi+i{y + he) - fi+i{y) for any y G (-6, -a - (/i^ H h he)) yields 

(A,) VKy) = (A/.)'7^+i(y + he) - (A,)'^/,+i(y) 

for any y € {—b, —a — {hd H h he)) and any k = — 1. The right-hand 

side is, however, nonnegative by (iv) of Lemma 2. This immediately yields 
that fe satisfies (A.l) for any k = 1, ...,£ — 1. 

Now, application of Lemma 2 to /2 implies that /2 is nondecreasing on 
{—b, —a — {hd + • • • + h2)). Because x, x + /ii G (—6, —a — hd — • • • — /i2) by 
assumption and 

Mx + hi) - Mx) = Ah,--- AhJix + hi) -Ah,--- AhJ{x) 

= Ah,---AhJ{x), 

it follows that A^, ■ ■ ■ Afi^f{x) > 0. Since A^^ - - - Ah^f{x) > can be shown 
along the same lines for any k = 1, . . . ,d— 1, the desired implication follows. 

□ 

Proof of (i) (in). First recall that if d = 1 and d = 2, respectively, 
d-monotone monotonicity of / reduces to / nonincreasing and nonincreasing 
and convex, respectively. It is therefore sufficient to restrict the discussion 
to the case d> 2. 

That (i) implies (ii) can be established by the midpoint theorem. Note 
first that if / is d-monotone on (a, b), then, for k = 1, . . . ,d — 2, f^^^ exists 
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on (—6, —a) and is nonnegative there. Consequently, for any k = 1, . . . ,d— 2 
and any x G (—6 — a) and any /i > so that x + kh ^ (—6, —a) there exists 
X* G + fe/i) so that {^hff{x) = P\x*) > 0. 

The reverse imphcation can be estabhshed using the same argument as 
Theorem 7, Chapter IV of [40]. To do so, pick an arbitrary x € (—6, —a). 
The assertion (iv) of Lemma 2 yields that, for any k = 1, . . . ,d — 1 and any 
h>0 such that x + kh e {-b, -a), (A,,)'=(/(x) - f{x - > whenever 
^ > is sufficiently small so that x — I G { — b, —a). By letting £ — > 0, it then 
follows that the left-hand derivative f'_ fulfills (A.l) for any k = 1, . . . ,d — I. 
The assertion (ii) of Lemma 2 in particular gives that f'_ is continuous on 
(—6, —a). Yet another application of Lemma 2, this time of the statement 
(iii) on /, then guarantees the existence of /' on (—6, —a). 

The same chain of arguments can be applied successively on /^'^^ In 
the last step, this yields that Z^*^"^) exists on (—6, —a) and fulfills (A.l) 
for k = 1,2. In particular, /("^"^^ is continuous, nondecreasing and con- 
vex. Put together, the derivatives f^^\ k = l,...,d — 2, exist on (—6, —a) 
and are continuous, nonnegative, nondecreasing and convex there. Because 
= the proof is complete. □ 

Remark A.l. After finishing this manuscript, Alfred Miiller (personal 
communication) brought the notion of d-convexity as defined and studied by 
[18, 31, 32] to our attention. It is not our aim to discuss the various concepts 
of higher order convexity here; nonetheless, we would like to mention the 
paper [3] containing in particular the following results: First, a continuous 
function / on (a, 6) which is weakly d-convex [meaning that {Ah)'^f{x) > 
for any /i > so that x + dh & (a, b)] fulfills A/^^ • • • A/j^/(x) > whenever 
x + J2i=ihi £ (o-jb). Second, if a function / is continuous and weakly d- 
convex on (o, 6), exists and is convex on (a, 6). These results would 

offer an alternative proof to parts of Proposition A.l. However, we favor our 
proof of the latter proposition, which is self-contained and more accessible. 
Furthermore, Proposition A.l relates d-monotonicity, which is a stronger 
concept than weak d-convexity, to nonnegativity of higher order differences. 
This link is absolutely central to the study of quasi-monotonicity of survival 
functions as delineated in the proofs of Proposition 2.2 and Theorem 2.2 
below. 

The fact that Proposition A.l characterizes d-monotonicity in terms of / 
rather than / may seem somewhat less elegant. The reason for this is that d- 
monotonicity is related to the existence of survival rather than distribution 
functions as described in Proposition 2.2, which we prove next. 

Proof of Proposition 2.2. Assume first that there exists a random 
vector X on MJ^ so that H{x.) = P(X > x) for any x € M!^. It is clear that 
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lim^^oo fix) = and /(O) = /(0+) = P(X > 0) € [0, 1]. Now, define f{x) by 
f{—x) for any x G (— oo,0] as in Proposition A.l and observe that for any 
X G ( — CX3, 0), any k = 1, . . . ,d and any /i > so that x + kh < 0, 

k 

^ "'"^ f{x + ih) 




X 

'I 



-X,<- + h 



n 



The last expression is, however, clearly nonnegative, in which case Propo- 
sition A.l implies that / is d-monotone on (0,oo). Because / is right- 
continuous at 0, it is d- monotone even on [0,oo). 

To establish the converse, note that / is in particular continuous on [0, oo) 
by Lemma 2. Therefore, H is right-continuous. Because condition (i) of 
Lemma 1 is clearly fulfilled as well, the only assertion which needs to be 
verified is the quasi- monotonicity of H{—x). Since H is right-continuous and 
places no mass outside of [0, oo)'^ and because l{x < 0} is quasi-monotone, 
this in turn amounts to the confirmation of Ah/( — ||x||i) > for any x € 
(— cxOjO)'^ and h > so that x + h < 0. However, this is immediately guar- 
anteed by (ii) of Proposition A.l upon setting x = ||x||i. Hence, the proof is 
complete. □ 



With Proposition 2.2, the proof of Theorem 2.2 is straightforward. 

Proof of Theorem 2.2. As already discussed in Section 3, C always 
satisfies the boundary conditions (i) and (ii) of Definition 2.1. Therefore, 
what remains to be shown is that C is quasi-monotone if and only if ip 
is d-monotone on (0,oo). However, because any Archimedean generator is 
continuous and fulfills ip{0) = 1, F given by F{x) = ^(x) for x > and by 
F{x) = 1 otherwise, is a continuous univariate survival function. 

Suppose first ip is d-monotone on (0,oo). Then, by Proposition 2.2, the 
function H defined as 

^(x)=V^(||max(x,0)||i) 

for any x € M'^ is a survival function in with continuous (identical) mar- 
gins F. Theorem 2.1 then ensures that the function C associated with H 
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via 

C{n) = HiF~\ui),...,F~\ua)) 

for any u G [0, l]'^ is a copula. Furthermore, F~^{ui) = mf{x ■.■ip{x) < u} = 
i/j~^{ui) for Ui € [0,1) and F~^{1) = mi{x : F{x) > 0} = — oo. However, 
max(F~^(l),0) = = and hence 

c(u) = ^{^-\ui) + ■■■ + V'-'M), u E [0, 1]'^. 

Conversely, assume that C is a copula. By virtue of Theorem 2.1 it there- 
fore follows that H{x) = C{F{xi),...,F{xd)), x G M"^, is a joint survival 
function with (identical) marginal survival functions F. Furthermore, it is 
easy to convince oneself that H['k) = ip{xi + • — \-Xd) for any x G M^. There- 
fore, H satisfies the hypothesis of Proposition 2.2 and it follows that ip is 
d-monotone on [0, oo). □ 

APPENDIX B: PROOFS FROM SECTION 3 

Proof of Proposition 3.1. Throughout, (y) indicates Stieltjes in- 
tegration; any other integral is to be understood in the usual Lebesgue- 
Stieltjes sense. The original result by [41] states that a function / is d- 
monotone on (0, oo) if and only if there exists a nonnegative, nondecreasing 
real function 7 on [0, 00) which is bounded below and such that 

Jo 

for any t G (0, 00). Furthermore, 7 is uniquely determined at its points of 
continuity and fulfills 7(0) = as well as 7(0+) = /(oo). 

Proof of (i). Suppose / is a Williamson d-transform of a nonnegative 
random variable X with distribution function F. Then /(0+) = /(O) = 1 — 
F{0) and lima;^oo f{x) = by dominated convergence. That / is d-monotone 
on (0, 00) is then clear from Williamson's result. 

To establish the reverse implication, observe first that in the situation of 
Proposition 3.1, Williamson's result ensures the existence of a unique right- 
continuous, nonnegative, nondecreasing function G which is bounded below 
and on M so that 

fit)= / (l-nt)t'dG(n) 
Jo 

for any t G (0, 00) and G(0— ) = G(0+) = as well as limu_^oo G{u) =p. The 
transformation formula for Lebesgue-Stieltjes integrals then yields 

f{t)= r(i--Y \v{u), 

Jo \ uj , 
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G(n) 



F{u) 



where i' is an image measure of the measure induced by G with respect to 
the mapping (j) given by (^(u) = l/u for u G (0, oo) and by otherwise. A 
possible distribution function of is 

0, if ti e (-00,0), 

_p-G{l/u-), iftiG[0,(X)). 

The disturbing fact that ly is not necessarily a probability measure can be 
easily remedied upon setting 

0, if -u € (-00,0), 

l-G(l/n-), iftiG[0,oo). 

It then still holds that f{t) = J^{1 - ^)+"^ dF{u) whenever t > but also 
that f{t) = 2IJdF(t) for any t G (0,oo)! Since both / and WaF are right- 
continuous at 0, the reversed implication is established. □ 

Proof of (ii). The fact that the relationship between / and F is one- 
to-one is clear from the original result by Williamson in combination with 
(i). Theorem 3 of [41], which is the inversion formula for Williamson d- 
transforms, further ensures that F is of the form (3.3) at all points of con- 
tinuity. Furthermore, the one-sided derivatives of (— l)'^~2y('^-2) exist ev- 
erywhere and are nondecreasing by Lemma 2. A version of the midpoint 
theorem for continuous functions with one-sided derivatives (see [25]) easily 

yields that is right-continuous on (0, oo) with left-hand limit . 

□ 
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Proof of Proposition 4.1. Throughout, let X ~ and U C be 
random vectors and be the Lebesgue measure on M'^. 

Proof of (i). It is sufficient to consider d = 2. In the case d>3 the situa- 
tion is simpler because then ip' exists everywhere in (0, oo) and one can argue 
by virtually the same arguments. Now, suppose H has a density. Then ob- 
viously P(X G (0,^-^(0))'') = 1. Further more, consider the transformation 
T:{0,i;-\0)y ^R"^ given by 

T(xi, ...,Xd) = mxi), . . .,ij{xd)), X G (0, V'-'(O))'^. 

Then T is injective and U = T(X). Because ip' exists a.e. in {0,ip~^{0)), T 
is a.e. regular (meaning that the set of points where T is not differentiable 
has Lebesgue measure zero). Furthermore, if ip'{t) exists, then necessarily 

To see this, assume for the moment the converse. Because ip is convex 
and decreasing on (0, V'~^(0)), V+ is nonpositive and nondecreasing. If now 
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< t < 'i/'~^(0) were such that i(^'{t) = 0, there would exist an e > so that 
^jJ'{x) = for almost all x G [t,ip~^{0) — e]. However, ip is absolutely contin- 
uous by assumption so il!'{x) = a.e. in [t,ip~^{0) — e] would imply that ijj 
is constant there, which is obviously not the case. 

Put together, D = {u e {0,lf:l\f^i = 0} is an empty set. 

In particular, Arf(L') = meaning that T is a.e. smooth. By means of Sard's 
theorem and the second transformation theorem (see [17], Sections 8.9 and 
8.10), it then follows that C has a density. 

The reverse statement can be established similarly by considering a trans- 
formation T{0, 1)'^ — > R'^ given by 

r(ni, ...,ua) = i^-Hui), . . . , V^-^M), X e (0, ly. 

Clearly, P(U S (0, 1)*^) = 1 and X = r(U) (for the latter claim, see the proof 
of Theorem 2.2). Again, one can easily convince oneself that T is a.e. regular. 
Furthermore, is (strictly) increasing, continuous and convex on (0, 1) 
and therefore in particular absolutely continuous in any [a,b] C (0,1). By 
the same argument as above, it can then be established that 'tp~^ (t) > 
whenever it exists. Consequently, T is a.e. smooth and X has a density 
again by Sard's theorem and the second transformation theorem. □ 

Proof of (ii). Suppose now ^p € "^d+i and let C be the {d+ l)-dimensional 
Archimedean copula with generator V'- Furthermore, fix a random vector X 
following the £i-norm symmetric distribution associated with C in the sense 
of Theorem 3.1. Arguments detailed in Section 5.2.3 in [8] then yield that 
all lower-dimensional margins of X have densities. In particular, this applies 

to Y = {Xi, . . . ,Xd)- However, Y = X and C has a density by means of (i). 
□ 

Proof of Proposition 4.2. Because of (i) of Proposition 4.1 and 
(iii) of Lemma 2, we only need to verify that the radial part of the ii- 
norm symmetric distribution associated with C has a density if and only if 
is absolutely continuous on (0,oo). Before we do so, observe that for 
X G [0, oo), Ffi may be rewritten as 

Fr{x) = Fq{x) {d-l)\ ' 

where Fq refers to the distribution function of the radial part of the ^i-norm 
distribution associated with the {d — l)-dimensional margin of C, that is, 

with {ui) H \- {ud-i)) ■ Furthermore, Proposition 4.1 guarantees 

that Fq is absolutely continuous. 

Assume first that Fr is absolutely continuous and consider [a, 6] C (0, oo). 

Because ft~^\x) = {-lf~^{d- lyS-^^^^^^^^^^ on [a, 6], /^"^ is absolutely 
continuous there. 
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To establish the reverse imphcation, observe that because /+~^ is contin- 
uous on (0,00), Fq is continuous on M and Fr is continuous at 0, then Fr 
is continuous in M. Now, fix [a,b] C (0, 00). By a similar argument as above 
it then clearly holds that Ff( is absolutely continuous in [a,b]. Because Fji 
is continuous and of bounded variation, it is absolutely continuous even in 
[0,6]. Since Fr is clearly absolutely continuous on (— c«,0], the absolute 
continuity of Fr on M is immediate. 

The last claim follows from the fact that if H is an absolutely continuous 
distribution function then its density h satisfies 

for almost all x € M'^. However clear this statement may seem, it is far from 
obvious; see [34], page 115 or [6]. □ 

Proof of Proposition 4.3. First recah that if U ~ C, then the ran- 
dom vector X = {tp~^{Ui), . . . ,tlj~^{Ud)) follows an £i-norm symmetric dis- 
tribution associated with C. Denote the radial part of X by and its 
distribution function by Ffj. Then, for any s G (0, 1], 

= p{J2X, = ^l^-Hs)^ =P{R = ij-\s)), 

where the last equality follows from (iii) of Proposition 3.2. Inversion (3.3) 
and the fact that ip^^^ is continuous for any k = 1, . . . ,d — 2 in turn imply 
that 

P(i? = ^-^(s)) 

= FR{^p-\s))-FR{iJ-\s)-) 

= {d-i)i (■v)-^+ (V' is)))- 

Regarding L{0), one similarly argues that P^{L{0)) = P{R > V^-^O)). Prom 
this it is first immediate that P'^(L(0)) = whenever ^~^(0) = 00. Suppose 
now ■i/'~^(0) < 00. Because all derivatives il>^^\x), k = 1, . . . ,d — 2, as well as 

■0+ ^\x) vanish for x £ [ip~^{0), 00), it ensues from (3.3) that -Fr(x) = 1 for 
any x S [ip'^ (0) , 00) . Therefore, 

p^{m)=p(j2i;-'m=i^-\o)y 
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Along the same lines as in the case of L{s), s S (0, 1], one reasons that 

P^(L(0)) = P{R = (0)) = ^"'^7i-"iy^^'"' ^-"'^ 
which concludes the proof. □ 

Proof of Proposition 4.6. First, recall that is Archimedean, as 
detailed in Example 2.2. Furthermore, ip^{x) = (1 — 2;)+ and Cl" coincides 
with the Frechet-Hoeffding lower bound. Therefore, assume that d > 3 and 
let be the generator of C. Because ~ (1 — 2;"^/'^'^""^)), the assertion 

^cL C can be shown by verifying that the function /(x) = 
is concave on (0, 00); see, for example, [27], Corollary 4.4.4. If ^p~^{0) < 00, 
it holds that f{x) = 1 for x G [^"^(0), 00) and it is sufficient to show that / 
is concave on [0,'ip~^{0)]. 

Assume first d > 4. In that case, / is twice differentiable on (0,^~^(0)) 
with 

/(2)(^) = __l_^(2)(^)^(^)-(<i-2)/(d-l) 

_ {d-2){4;'{x)f -{d-l)^P^^^x)i;{x) 
(d-l)2V;(x)(^-2)/('^-i)+i 

To establish that f^^x) < for x G {0,ip~^{0)), we consider the matrix 

^l^{x) 

A = 




V d-1^'^''^ (d-l)(d-2) 

and show that it is positive semi-definite for all x G (0, ?/'~^(0)). Since "0 
is d-monotone, the dominated convergence theorem ensures that, for x G 

(o,rno)), 

roo 1 / ™\ d—l — k 

4,i^){x) = {d-l)-..{d-k)j^ dFnit), 

k = l,...,d-2, 

where is the distribution function of a nonnegative random variable 
whose Williamson d-transform is ip. Now, fix Oj G M, i = 1, 2. Then 



^ ^ /'OO / \ ct~ i 



i=l i=l 



38 



A. J. MCNEIL AND J. NESLEHOVA 



+ 2aia2 

+ ia2? 
1 - 



oo ]^ 

t 

oo ^ 
^ 

1 



1 



d-2 
+ 

d~3 



dFnit) 



t 

ai(l 



dFnit) 



+ 



+ 



02 



dFij(t) > 0. 



Consequently, |^| > 0, which in turn imphes f^'^\x) < for x G (0,'i/'~^(0)). 

For d = 3 we can use the fact that Tp^^^ exists everywhere in (O,'0~"^(O)). 
Consequently, f^\x) exists for any x S (0,-0""^ (0)) and 

(2) (d-2)(V.^(x))2-(d-l)^f (x)^(x) 

■1+ W (^_i)2^(3.)(d-2)/(d-l)+l 

Concavity of / then follows from < on (0,'i/'~^(0)); see Theorem 1 of 
[25]. Now consider h>0 and x G (0,-0"^ (0)). Then 

Tp'i^x + h) — iIj'{x) 



&(1A)(1 -{x + h)/t) dFnit) - /r (1A)(1 - ^/t) dFnit) 



h 



f dFnit) + f -(l-^) dFnit) 

J(x+h,oo) t^ J{x,x+h] th\ t J 

Because ^(1 — |) < ^^^.''^^■^ for t € + /i], the dominated convergence 
theorem yields 



V.f(x)=lim^:^^^±^l^^ = 2 



^ dFnit). 

{x,oo) I 



The inequality f^\x) < for x G (0,V'"^(0)) can now be established by 
exactly the same arguments as in the case d > 4. □ 
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